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一 类 求解 带 随 机 成 本 的 生产 运输 问题 的 线性 通 近 方法 * 


万 中 ， 江 P, HÆR, TEH 
(中 南大 学 数学 科学 与 计算 技术 学 院 ， 长 沙 410083) 

摘 要 : 生产 和 运输 集成 计划 问题 在 许多 工业 工程 领域 都 普遍 存在 。 要 给 出 最 优 的 生产 和 运输 计划 就 必须 
考虑 实际 工业 管理 过 程 中 存在 的 不 确定 性 因素 。 本 文 研究 了 生产 厂家 的 生产 能 力 、 商 家 的 需求 量 
和 单位 运输 成 本 等 因素 为 随机 变量 情况 下 的 产品 生产 与 运输 成 本 问题 ， 建 立 了 该 类 问题 的 随机 优 
化 模型 。 在 一 定 的 假设 条 件 下 ， 推 导 了 所 建 模型 的 确定 等 价 类 。 基 于 问题 的 结构 特征 ， 提 出 了 求 
解 生产 和 运输 计划 的 一 种 线性 逼近 方法 ， 数 值 例子 表明 该 种 方法 的 应 用 前 景 。 
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其 目标 就 是 要 极 小 化 生产 和 运输 成 本 ， 其 中 目标 函数 中 的 第 一 个 和 式 (线性 部 分 ) 表示 总 的 运输 
成 本 ， 第 二 个 和 式 (和 性 成 本 ， 可 分 离线 性 或 非 线 性 部 分 ) 表示 生产 厂家 的 总 生产 成 本 。 其 它 
记号 的 意义 如 下 ，!: 所 需要 运输 的 产品 种 类 ; m: 生产 厂家 的 个 数 ，n: 商家 的 个 数 ， wki: 
第 i 个 生产 厂家 生产 产品 上 的 产量 ，gki(wki): 第 i 个 生产 厂家 生产 wki 所 需 的 生产 成 本 ;cpij: 
生产 厂家 i 送 到 商家 j 的 产品 上 的 单位 运输 成 本 ; mug: 生产 厂家 i 给 商家 j 提供 产 品 的 送 货 
量 ， 也 是 决策 变量 。 

者 用 ou 表示 第 i 个 生产 厂家 生产 或 供应 产品 上 的 能 力 ，i = 2 2 …… m: 0x; 表示 第 了 个 商 
家 对 产品 上 的 需求 量 ，j = 1,2,… ,n， 且 设 api 和 bx; 是 都 服从 正 态 分 布 的 独立 随机 变量 ， 即 


aki ~ N (Hki, o2), bk; ~ N (uk, vi), 1—1,2,.-.,m J212,.-.,n k=1,2,... ,1, 
则 约束 条 件 


Ds Do 


-1j-1 
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可 转化 成 上 述 模 型 的 第 一 个 约束 ， 其 中 ax; 是 此 随机 不 等 式 得 到 满足 的 置信 水 平 。 同 理 约束 条 
件 


m 


1 
Tkij 之 5 bki 
k=1 


k=1 i=1 


可 转化 成 上 述 模型 的 第 二 个 约束 ， 而 B; 是 相应 的 置信 水 平 ，g-!(.) 是 标准 正 态 分 布 的 分 布 函 
NUN SCA. DRE, aki, Bkj € [0, 1]. 

ZXR T FÆ ari 和 bkj 看 作 随 机 参数 外 ， 考 虑 到 实际 中 产品 类 型 ， 油 价 涨 中 ， 公 路 交通 状况 
以 及 司机 技能 等 各 种 因素 均 会 影响 产品 的 单位 运输 成 本 cki;}， 因 而 视 随 机 参数 处 理 更 合适 ， 我 
们 将 在 文献 [1-3] 的 基础 上 重新 推导 这 类 随机 模型 的 确定 型 等 价 类 (EDF)， 并 设计 有 效 的 求解 
算法 。 


2 多 产品 运输 问题 的 CCP 模型 


在 模型 (1) rn, did 


D ete o Cerati ela Gimis e a Cmn yr i ,Comny ee s Elan s 
X? = add oet dic Flacc odictor dito dro ludi t llam e aedi: 
1 1 0 0 0 0 
0 0 1 1 0 0 
U = 0 .0 0 .0 .1 .1 |, EE=(,1,. ,1)ixdm), 
—1 0 -1 0 —1 0 
d oom EE DIU me st aT 
人 
D' = (pa 一 aag (01), 777, fim — Tim (03m), Ha 一 aag (o), o, 
Him — Omp (Om), >, -vii — e (Bri) c, Vin — m Ó (Bin) 


— ui — 019^! (E), +, -Vin — emÓ (Bin)), 


HH U fü D 2 IDE (lm + In) x (Imn) 和 (lm + Un) x 工 维和 矩阵 。 此 时 模型 (1) 可 写 出 如 下 向 量 形 
式 
min F(X,G)-C'X + EG 
st UXXD, (2) 
X 20. 
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利用 机 会 约束 规划 方法 5; 引 ， 模 型 (2) 又 可 以 转化 成 下 列 形式 
min € 
st. P(F(X,G) <&€)=7, 
UX <D, 
X 20. 


设 C 的 分 布 函数 已 知 ， 且 数学 期 望 为 C， 协 方差 矩阵 为 Vmn) x mn) y € [0,1] 是 给 定 置 信 水 
^F. HU (3) 对 应 的 确定 型 等 价 类 为 


(3) 


min F(X,G)-ÓC'X +4! (y X'VX + EG 
st. UX < 也， (4) 
X. 


下 一 节 我 们 将 设计 求解 (4) HAARA, REEN. 


3 ”最 优 解 的 性 质 和 线性 逼近 方法 


模型 (4) 求解 的 主要 困难 在 于 单位 运输 成 本 C 是 随机 参量 。 即 使 知道 C 所 服从 的 分 布 ， 目 
RRHH o VXV X 和 EG 都 是 非 线性 的 。 为 设计 求解 该 问题 的 有 效 算法 ， 我 们 先 研 究 
这 类 问题 最 优 解 的 性 质 。 

引 理 1 s, sz 是 任意 两 个 不 同 的 随机 变量 cl, cz 的 方差 ，s12 是 它们 的 协 方差 ， 则 有 


11a 1c 
—$154 < $12 < Si S2. (5) 


证 明 对 任意 非 零 二 维 向 量 w， 因 为 


( ) E(cı Ec)? 五 (ci s Eci)(c2 B Eca) Ql 
pa 五 (cl x Eci)(co — Ec») E(co -— Ec2)? 


2 
ES E( 》 "aici(c 一 Ba) >0, 
i=1 


所 以 和 矩阵 
E(ci — Ea)? E(ci — Eci)(c2 — Ec) 
E(ci — Eci)(co — Ec2) E(co — Ec2)? 
是 正定 矩阵 ， 其 所 有 主子 式 大 于 零 ， 从 而 其 对 应 的 行列 式 大 于 零 ， 即 sls? > VSi12。 两 边 平方 
可 得 到 不 等 式 (5)。 


引 理 2 对 任意 非 负 实 数 zl £2, DA 有 如 下 不 等 式 
(ri-- 22. +n)? <L? +r +r (6) 


成 立 。 
根据 上 述 两 个 引 理 ， 可 以 证 明 如 下 结论 。 
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引 理 3 设 卫 是 问题 (4) 中 的 决策 变量 ， 则 


imn 
C'X —ó (y X sêt + EG < C'X có (yv X'VX + EG 

g—1 

imn : 
< CX (y) 》 sin, + EG. (7) 
q-1 
id 
imn imn 


F(X,G)= C'X 4 $71(4) X six, -HEG, FYX,G)-C'X—4 (y) 5 siz, + EG. 
gq=1 q=1 


E S 是 模型 (4) RRE, frRI Sy 分别 是 问题 


min fı(X,G) 
s.t. UX <D, (8) 
X>0 
和 问题 
min f(X, G) 
st UX <D, (9) 
X20 
的 最 优 解 ， 则 根据 引 理 3 有 
ISFE (10) 


假设 1 问题 (8) 存在 一 个 有 限 最 优 解 X?。 
定理 1 如果 假设 1 成 立 ， 则 问题 (4) 的 最 优 解 也 是 有 限 的 。 
证 明 假设 问题 (4) 的 一 个 可 行 解 Xo 满足 


imn 


C Xo $^ (y) VXV Xo + EG 2 O Xt 4 $7 (y) Y s?a} + EG. (11) 
q—1 
因为 
imn 1 
$7 (y) / X)V Xo + EG < 9 (y) Y sia, + EG, (12) 
q-1 
所 以 
RI X? < C'Xo- à! (q)/ XV Xo + EG € R} Xo. (13) 
上 式 中 
imn imn 


RIX} =Ö'X} +H) Y sbX; + EG, RQXo-Ó Xo 47 (v) Y. s? Xo, + EG. 
q=i 


q=1 
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由 此 可 得 R X? < RiXo。 这 与 假设 1 矛盾 。 因 此 对 任意 一 个 可 行 解 Xo， 下 列 不 等 式 成 立 
imn 


C Xo 47 (y) XIV Xo + EG < O' Xt 4- 97! (y) V s? Xt, + EG. (14) 


q=1 


假设 2 问题 (9) 存在 一 个 有 限 最 优 解 X3。 

定理 2 如 果 假设 2 成 立 ， 则 问题 (4) 的 最 优 解 也 是 有 限 的 。 

证 明 过 程 类 似 定理 1 的 证 明 。 

根据 上 述 讨 论 ， 我 们 构造 求解 (4) 的 如 下 双 目 标 线性 规划 子 问题 


= imn 1 lmn 
min fı(X)=C'X +y) 35 siz,4 EG = 》 Cito 
q=1 q=0 


pe lmn 1 lmn 
min fo(X)= C'X—4 !(y) 35 sáz4 + EG = Y, eaqtq (15) 
9 一 1 q=0 
s.t. UX <D, 
X 20. 
通过 求解 子 问题 (15)， 我 们 将 得 到 问题 (4) 的 最 优 解 的 近似 解 回 。 因 为 双 目 标 问题 (15) 的 每 个 
次 优 解 都 可 能 是 模型 (4) 的 一 个 最 优 解 ， 下 面 研 究 它们 之 间 的 关系 。 
定理 3 WR X* 是 模型 (15) 的 一 个 次 优 解 ， 且 存在 一 个 w*， 满 足 


f(z") = w* fi(z*) + (1 ~ w*)fa(z") = F(z"), (16) 


则 0.5 < w* < 1. 
证 明 由 (16)， 我 们 可 以 得 到 


w* = (C'X* ce! (y)vxnVX* — fa(X*)/(RQC) — RX). (7) 
$s? = (s3, s8), 于 是 有 
Rh(z*) 2 C X* + 97! (g)s? X* + EG, (18) 
fa(7”) = C'X* — $71 (yet X* + EG, (19) 
所 以 
w* = (s2 X* + (X"V X")3)/(252 X"). (20) 


Dy X"VX* < s3X*， 所 以 w = X"VX*/253 X* < 0.5。 于 是 w* = cuX TF 05fB i 
过 1。 定 理 3 得 证 。 
定理 4 模型 (4) 的 最 优 解 属 于 模型 (15) 的 次 优 解 集 。 
证 明 设 模 型 (4) 的 目标 函数 的 最 优 解 是 X*， 且 存在 一 个 w*， 使 得 X* 在 等 式 (7) 成立， 
又 设 存在 X' 是 下 列 模型 的 解 
min fi(X)+fo(X) = 20'X+2EG 
st. UX <D, 
fi(z) € fi(z"), (21) 
falx) € fa(x*), 
X 20. 
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在 模型 (21) F, fi (x) 和 户 (z*) 是 模型 (4) 在 最 优 解 处 的 函数 值 ， 而 模型 (21) 的 最 优 解 六 
是 属于 模型 (15) 的 有 效 解 的 解 集 之 中 。 因 为 XX' 是 模型 (4) 的 可 行 解 ， 且 满足 下 列 两 个 不 等 式 


filz’) < fi(z’), (22) 
fala’) € fia"). (23) 
* 
F(z) = wi filz) + (0 — wi) f(z), (24) 
其 中 
Wi 二 一 xvi wo 一 一 w 
1 二 了 (253 X) s 2 1; 
因此 可 得 
F[fi(z^), fa(2^)] > F[fi(z^), fa(2*)]. (25) 
于 是 
F(X^) = f(X’) 2 F(z*) = f(z*), (26) 


又 因为 z* 是 模型 (4) 的 最 优 解 ， 所 以 z* = X'。 于 是 定理 4 成 立 。 
314 BERR% f(x) EIT H = {UX < D,X>0}+ 上 的 最 优 解 不 小 于 户 (z), 户 (z) 在 
可 行 域 瑟 上 的 算术 平均 数 ， 即 


FE) 2 51 2) 2 5 (G7) + i). 


证 明 从 略 。 

定理 5 在 可 行 域 忆 上， 相对 于 户 (z) 而 言 ，f(z) 的 最 优 解 更 加 接近 fi (v) 的 最 优 解 。 

证 明 令 Di -—|f(z")- ffl, Do —|f(z*) - f: RGEUEBH Di € Do 成 立 。 假 设 Di € Ds 不 
R, M Di > Do. MAJE) - f) > lz l 因为 fo(z*) < «fm. 
ft» f(z*), BIB —f(z*) - ft > f(z*) - f. XXE 2f(z*) < +f Awa F F 
是 假设 不 成 立 ， 即 定理 5 成 立 。 

定理 4 表明 模型 (4) 的 解 可 以 通过 求解 模型 (15) 的 次 优 解 集 求 出 ， 但 求解 模型 (15) 的 次 优 解 
集 需 要 花费 一 定 的 时 间 。 我 们 可 以 利用 一 个 折 中 函数 来 说 明 求解 的 过 程 在 何 种 状态 之 下 可 以 停 
止 ， 这 个 函数 可 以 定义 为 

imn imn 
g(z) = tUbnm ^ b :( Yours f2 ) Inin -( i A Js (27) 
> 2 P qn (X 2 P qct 


wi, w 是 用 (x) 和 f(z) 规范 化 后 的 权重 系数 ， 且 ui + ws = 1. H (27) R4 g(z) 仍然 是 一 个 线 
性 函数 。 若 将 g(z) = 0 为 约束 条 件 (不 妨 称 之 为 折 中 条 件 )， 则 最 优 解 能 够 保证 g(z) = 0, HE 
这 种 情况 下 在 寻求 模型 (4) 的 次 优 解 的 过 程 当 中 折 中 条 件 的 个 数 是 有 限 的 。 


4 ”算法 描述 和 算 例 


根据 上 述 讨 论 ， 我 们 可 以 设计 如 下 求解 原 问题 的 算法 。 该 算法 由 两 部 分 组 成 。 
第 一 部 分 如 下 : 
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步骤 1 写 出 线性 规划 模型 (8), (9) 的 具体 表达 式 ; 

步骤 2 求解 线性 规划 模型 (8), (9) 的 最 优 解 m1, c5. 

第 二 部 分 如 下 : 

YE ui, wu 分别 是 权重 系数 书 的 一 个 下 界 和 一 个 上 界 ，e 是 精确 度 ，k 是 确认 模型 (4) 的 次 优 
解 的 计算 次 数 。 依 次 执行 如 下 步骤 : 

步骤 1 Bk = 0; w= 0.05, wu = 0.99, w = wu, Dw = 949, e = 0.00001 和 及 = f(x?)， 然 
后 执行 步骤 ?2; 

步骤 2 设 大 = 大 十 1， 写 出 g(z) 的 具体 表达 式 ， 且 9(z) = 0, 


imn imn 


(z) = f MÀ -fi 
B (Ese! rl SS | ) GPS ( » ciata : ) 
q—0 
再 执行 步骤 3; 
步骤 3 求解 如 下 问题 
min fi(X) 


(28) 
st. UX «D, g(z)-0, X20. 


H X} 表示 模型 (28) 的 最 优 解 ， 令 及 = f(X), WITHER 4: 

步骤 4 WR Fe > Fk- Ruwi = wh Wu = Wha Dw = iDw,wi = wa — Dw, WiTH 
骤 2; 否则 执行 步骤 5; 

步骤 5 WR Fk- 及 -<e， 则 达到 近似 最 优 解 ， 否 则 执行 步骤 6， 

步骤 6 计算 ww = wi 一 Dw, wh = w。 如 果 w < 0.5， 令 Du = 3Dw， 执 行 步骤 1; 否则 执 
行 步骤 2。 

利用 上 述 算 法 我 们 求解 一 个 含有 4 个 工厂 ，3 个 销售 地 ，2 种 产品 ， 且 生产 力 ， 需 求 量 和 单 
位 运输 成 本 都 是 随机 的 生产 和 运输 计划 问题 ， 其 中 取 ~y = 0.95， 两 种 产品 的 单位 运输 成 本 期 望 
和 方差 如 下 


Č = (22, 24, 31, 33, 35, 28, 35, 40, 37, 27, 39, 31, 37, 31, 35, 31, 43, 28, 35, 27, 26, 29, 35, 35), 
S = (1,2,2,1,2,1,1,1,2,1,1,1,1,2,1,2,1,1,1,1,2,2,1,2). 


车 令 思 性 生产 成 本 函数 gp; = 200 x (awki)06， 利 用 Lingo8.0 求 解 可 以 求 出 模型 (8), (9) 的 
fff = 29494.8, fj = 27763.7。 再 利用 第 二 部 分 的 算法 分 别 取 不 同 的 权重 系数 可 计算 
出 产 =29013.5， 当 7 取 不 同 的 值 时 ， 求 解 的 结果 如 表 1。 


表 1: 不 同 置信 水 平 下 线性 逼近 法 的 数值 结果 


y $^() f 
0.96 1.755 29368.5 
0.97 1.885 297721 
0.98 — 2.055 29813.5 
30186.4 


y $ (9 f y $70) J 
0.50 0 26150.2 0.90 1.285 27709.7 
0.80 0.845 26359.4 0.92 1.405 27990.4 
0.85 1.035 27247.7 0.94 1.555 28128.3 
0.88 1.175 27377.9 0.95 1.645 29013.5 
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Linear Approximation Algorithm for Production and Transportation 
Schedule Problems with Stochastic Cost 
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Abstract: The integrated production and transportation scheduling problems are very common in 
many fields of industrial engineering. An optimal scheme of production and transportation should be 
involved with some nondeterministic behaviors existing in the industrial management. In this paper, 
we construct a stochastic optimization model for production and transportation scheduling problem, 
where all of the capacities of production and transportation, the cost and the demand are stochastic. 
Under some settings, we derive the deterministic equivalent formulation of the original problem. Based 
on the structural features of the model, we have developed a linearization approximation algorithm 
for this class of problems. Numerical examples show that both the proposed model and the solution 
method are promising. 

Keywords: production and transportation scheduling problem; stochastic optimization model; lin- 
earization approximation method; deterministic equivalent formulation 
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